Transfer and storage of vortex states in light and matter waves. 
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We theoretically explore the transfer of vortex states between atomic Bose-Einstein condensates 
and optical pulses using ultra-slow and stopped light techniques. We find shining a coupling laser 
on a rotating two-component ground state condensate with a vortex lattice generates a probe laser 
field with optical vortices. We also find that optical vortex states can be robustly stored in the 
atomic superfluids for long times (> 100 ms). 
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One of the most dramatic developments in recent ex- 
periments on atomic Bose-Einstein condensates (BECs) 
is the creation and study of quantized vortices and vortex 
lattices 0,0- In this Letter we show that light modes 
with orbital angular momentum can be efficiently stored 
in superfluid vortex configurations and later rewritten 
back onto light fields. This also provides a feasible mech- 
anism to transfer BEC vortex lattice ground states to 
light pulses. The strong light-matter coupling can be im- 
plemented using the recent technology of ultra-slow and 
stopped light pulses in atomic vapors based on 

the method of electromagnetically-induced transparency 
(EIT) Q. In these experiments, an optical field's ampli- 
tude and phase pattern was written onto superpositions 
of atomic states, stored for several milliseconds, and then 
re- written back onto the original light field and output. 

One of the unique characteristics of slow light is the 
strong yet coherent action of the matter fields on the 
light fields We here study the robust storage and 
transfer of coherent information in terms of vortex config- 
urations, for experimentally feasible parameters of a two- 
component 87 Rb BEC. We find that shining a coupling 
laser pulse can generate a probe pulse which contains op- 
tical vortices corresponding to the vortices originally in 
the atomic fields. If vortices are present in both BECs, 
the vortex circulation of one BEC component is reversed 
in the probe field. This light field could potentially be in- 
put into a second BEC which is optically connected with 
the first BEC (even if it is spatially distant), allowing 
transfer of vortex states between BECs. 

The strong coupling between the quantized optical and 
BEC vortices is an interesting phenomenon in its own 
right, as a demonstration of nonlinear superfluid- light op- 
tics and may be useful in imaging vortex lattices. But it 
also points to possibilities for information processing and 
storage. Experiments have seen advances in the genera- 
tion of Laguerre-Gaussian (LG) modes of light with or- 
bital angular momentum as well as applications of these 
modes to quantum information It may be possible 
to utilize modes with several different angular momen- 
tum quantum numbers to build a quantum information 
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FIG. 1: Writing an atomic vortex lattice onto a light 
field. The optical density (top row) Di = Na f dz\i/ji\ 2 
(black=0, white=55) and the phase 4>i in the z — plane 
(bottom row) for a vortex lattice (a) in |1) and (b) in |2). 
(c) n p (the number of probe photons output per cross sec- 
tional area a) upon a switch-on of the coupling field (black=0, 
white=25) and the phase profile of 4> p , which contains vortices 
of both circulation, (d) Same quantities as (c) for a non- 
rotating ground state BEC in |1), resulting in vortices of only 
one circulation, (e) A scatter plot of the transfer efficiency 
T — n p /u2, where ri2 is the number of \2) atoms per a before 
the switch-on, versus /13-D1 for the case in (c) (open circles, 
red online) and for (d) (closed circles, blue online). The solid 
curve shows the estimate 1 — (1 + /i3-Di)~ 1,/2 . 



architecture based on a larger alphabet than the tradi- 
tional two-state systems [lOj. However, optically based 
quantum information schemes suffer from the difficulty 
of trapping and storing optical fields. 

Here we perform a comprehensive analysis of how LG 
beams can be written into BECs, stored in the atomic 
fields, and later rewritten onto light fields. We find one 
can choose parameters so the coherent BEC dynamics do 
not dissipate the information during the storage time, al- 
lowing storage fidelities on the order 70% for several hun- 
dred milliseconds with typical parameters. This fidelity 
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can be further improved by using BECs with larger opti- 
cal densities. Moreover, the phase information in vortex 
lines can exhibit very long life times of tens of seconds 
or can even be energetically stable 0. While similar 
LG beams have been proposed as a method to generate 
vortices in BECs |ll| an d m experiments vortices have 
been imprinted on BECs 2] , our scheme here is very dif- 
ferent: In the nonlinear light-matter coupling both am- 
plitude and phase are robustly and simultaneously ex- 
changed, allowing the controlled storage of information. 
Moreover, slow light and EIT are now inspiring impor- 
tant coherent matter wave technology applications, e.g ., 
BEC shock waves IgJ, p robing decoherence in BECs [12j , 
optical black holes Il3 l , storage and production of non- 
classical states 0, Il5j |. quantum information processing 
flif , lo w lig ht level nonlinear optics ^3 > an( i optical pro- 
cessing pL8| . 

We consider two BECs of Rb-87 atoms, in stable 
internal states |1) = 1551/2, F = 2,M F = +1) and 
|2) = 1551/2, F — l,Mp = —1), with macroscopic wave- 
functions tpi and ifa. These are connected, by reso- 
nant, +z propagating coupling and probe light fields, Q c 
and £l p (each of wavelength A = 780 nm), to the ex- 
cited state |3) = |5Pi/ 2) -F = 2,M F = 0), which de- 
cays at T = (2tt) 6 MHz, forming a A three-level struc- 
ture. Here the Rabi frequencies fl{(r, t) = ~di 3 • Ej(r, t) 
(p = 1, c = 2) arc defined in terms of the atomic dipole 
matrix elements and the slowly varying envelope 
(SVE) of the electric fields Ej (with the rapid phase rota- 
tion at the optical frequencies and optical wavenumbers 
factored out). The light propagation then follows from 
SVE approximation to the Maxwell's equations [19j: 

( d 1 d\ ik . _ 

Here the SVE of the polarization is ■ Pi 3 ~ 
f i3 V 2 Nilj*ip 3 , N is the initial total number of BEC 
atoms, /13 = 1/4 and ^23 = 1/12 are dimensionless os- 
cillator strengths, V — \J heoaT/k is the reduced dipole 
matrix element and a = 3\ 2 /2ir is the resonant cross- 
section. The BEC wavefunctions ipi , "02 evolve according 
to generalized Gross-Pitaevskii equations (GPEs): 

ihi>i = {H Q +J2 Uij m 2 )A + m*^ , (2) 
3=1,2 

where Hq = — h 2 W 2 /2m+V, with a perfectly overlapping 
trapping potential for both states V(r) = y [ui 2 (x 2 +y 2 ) + 
<jJ 2 z 2 ]. In the interaction coefficients Uij = AttNH 2 aij/m, 
the a,ij are the scattering lengths between atoms in states 
\i) and \j). For 87 Rb, 012 = 5.5 nm and an : a.12 : 022 - 
0.97 : 1.0 : 1.03 [2(|. Often in practice, and throughout 
this Letter, we will be working in a regime where, when- 
ever the light fields are present, the atomic internal state 
dynamics and light field couplings are on a much faster 
time scale and dominate the external dynamics. In the 



context of slow light experiments we can adiabatically 
eliminate ip 3 ~ -i(fl p ipi + n c ip 2 )/T in Eqs. (|Tl2"j> [lif . 
Then the last term in Eq. results in both coherent 
exchange between |1), |2) as well as absorption into |3). 
In our model, atoms which populate |3) and then spon- 
taneously emit are assumed to be lost from the BECs. 

Vortex lattices can be produced by imparting angular 
momentum to the system with a time-dependent rotat- 
ing potential. We numerically find the ground state of 
a two-component BEC rotating along the z axis at the 
speed Q — 0.3 uj r by evolving the GPEs (0 in imagi- 
nary time in the absence of the light fields in the rotat- 
ing frame, obtained by replacing Hq by Hq — u)L z , with 
L z = ih(yd x — xdy). At every time step, we separately 
normalized the wavefunctions to fix the atom number in 
each component. The full 3D integration, without im- 
posing any symmetry on the solution, is performed in a 
'pancake-shaped' trap uj r — 0.1w z = (2tt) 10 Hz, with 
N = 2.76 x 10 5 , on a spatial grid of 256 2 x 32. 

In Fig. H w e present a calculation of such a two-species 
lattice. The phases of the wavefunctions indicate a lat- 
tice of singly-quantized vortices in each component, with 
the positions of the vortex cores in the two species offset 
from each other. The filling of the vortex cores by the 
other BEC significantly increases the core size, as com- 
pared to vortices in a single-component BEC, and makes 
them observable even without a ballistic expansion. The 
3D results are qualitatively similar to the previous 2D 
calculations [2l]]: The rectangular vortex lattice pattern 
is recognizable close to the trap center, but becomes dis- 
torted close to the BEC boundaries. The total density in 
the center is 6.9 x 10 13 cm~ 3 . The vortex core positions 
(and (f>i,(f>2) change very little along z. 

To model the writing of this vortex lattice to the probe 
field flp we numerically solved Eqs. 1)112(1 when a coupling 

field with a peak input value of f2c = (27r) 4 MHz was 
switched on suddenly (in a time 0.8 (is). When this hap- 
pens the coherence V'iV'2 acts to generate a probe field 
Qp in such a way that the contributions to the polariza- 
tion in Eq. cancel out, and a dark state is formed. 
The resulting output intensity and phase pattern of fl p 
is shown in Fig. ^c). 

To understand the results, we note that solving 
Eq. ftp, ignoring the negligble z variation of ^1/^2, 
shows the fields eventually reach the asymptotic values 

n p = -4 in) (/i3^ 2 )/(/i3|^i| 2 + /23I02I 2 ) and n c = 

Q p ipi/ip2- Note that f2 p — > whenever either ~ * 5 
giving vanishing intensity at all the phase singularities. 
The phase 4> p is determined by the relative BEC phase 
4>p = 02 — 4>i + n (the coupling field acquires no phase 
shift here) . Thus 4> p contains singularities with one circu- 
lation at the location of vortices in ip2 and of the opposite 
circulation at vortices in 

In the columns where l^/V'il ^ 1 the coupling field 
coherently generates the probe field via coherent trans- 
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fers from 1 2) and 1 1) so that the number of probe photons 
output per area a, n p , equals the original number of |2) 
atoms per a, n 2 . Deviations from this ideal limit are due 
to absorptions in a 'preparation' region (equal to one op- 
tical depth) before the dark state is established. When 
the optical depth of |1), /13-D1 » 1, this loss is relatively 
small. On the other hand, near the vortex cores of |1) 
where tpi — > there is not enough coherence to gener- 
ate ri2 photons. Instead f2 c is attenuated by absorption 
and eventually the dark state is established primarily via 
incoherent spontaneous emission events. We found that 
the transfer efficiency T = n p /u2 for each column can be 
determined by the optical depth, agreeing very well with 
the prediction 1 — l/V/13^1 which approaches unity as 
Di — ► 00; see Fig. [He). The overall output efficiency 
integrated over x, y is 35%. 

Motivated by the fact that the fidelity improves with 
D\ we also considered coupling the ■02 vortex lattice of 
Fig.[IJb) to a vortex- free, non-rotating ground state BEC 
in |1), with N = 3 x 10 6 [Fig.Hfd)]. There the peak opti- 
cal density /13-D1 is 70 and is > 35 throughout the region 
occupied by |2). To a good approximation, the probe 
intensity simply reflects ip% [see Fig. ^b)] and only con- 
tains vortices of one helicity. The fidelity is now higher 
[Fig. Die)] with the overall efficiency 85%. 

Having determined that vortex states can indeed be 
written from atoms to light fields we now consider the 
possibility of storing information, input in the form of 
vortices in light fields, in a two-species BEC. For this 
study, we switch the roles of the two stable 87 Rb states 
so |3) is coupled to |2) = |F = 2,M F = +1) by fl c 
and to |1) = 1 1 , — 1) by Q p . A BEC is initially in the 

(non-rotating) ground state of 1 1) , labelled ip^ (r), in an 
isotropic trap uj z = ui r = (27r)21 Hz, with N = 4 x 10 6 , 
Rtf — 25 /j,m, peak density of 1.5 x 10 14 cm' 3 , and a 
peak optical density of /13D1 = 213. The initial coupling 
field is at «i in) = (2tt) 8 MHz [see Fig. 0a)], while the 
LG probe input reads (in cylindrical coordinates 23] ): 
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FIG. 2: Inputting, storing, and outputting an optical 
vortex in a BEC. (a) A spherical BEC (grey sphere) is 
originally illuminated by a c.w. coupling field £l c . A probe 
field pulse in an LG ('doughnut') mode (black torus), is then 
input and stopped when it has propagated to the center, (b) 
The effective potential seen by tp2, V"2_ c ff(r) in units of tkJ z , 
during the storage for m = 1, with a minimum near r — 6fim, 
z = 0. (c) The number of atoms in |2), N2, relative to the 
number of input photons Np lu ^ (solid curve) and the number 
of atoms lost due to spontaneous emission (dashed curve) dur- 
ing the probe pulse input and coupling switch-off (the vertical 
dotted lines shown the points at which the coupling field is 
90% and 10% of its full value). These are constant during the 
storage time (the time break). After the time break we then 
plot the same quantities during the coupling switch-on and 
probe pulse output, (d) The relative density profile \ip2\ 2 /p 
during the storage at the times indicated. The dotted curve 
indicates the BEC boundary. The dotted lines shows the cuts 
plotted in (d). (e) \ip2\ 2 /p and phase of cf>2 at the same lines 
along the cuts indicated in (c). (f) The corresponding output 
fidelity T upon switch-ons after various storage times r s t, for 
m — 1, 2, (solid, dot-dashed, and dotted lines, respectively). 



to absorption, and so the number of spontaneous emis- 
sion events Moss reaches about 0.2lN p n ^ at this point. 

This pulse is sufficiently weak |f2 p | <C |f2c n ^| that the 
coupling field il c is nearly homogenous across the BEC 
3.83 x 10 4 photons and propagates slowly, with and tpi ~ ?p[ G ^ during the propagation. Then -0 2 is co- 
herently driven into the dark state ip2 — — (Cl p /£lc ) 
and so acquires the amplitude and vortex phase pattern 
of the input LG mode. At t = 0.6 /xs, f2 c is switched off 
in about 0.2 yus, and N2 and A^ oss are unaffected by this 
rapid switch-off; see Fig.^c). Both fl c and fl p smoothly 
ramp to zero intensity and the LG mode is then stored 
in ip2- Figures EI d-e) show the relative density profiles 



w 



|m| 



(3) 



where w is the beam waist size. Each probe photon in 
the LG mode with |m| > 1 contains a vortex at r = 
with a vanishing intensity and m units of orbital angular 
momentum. 

We consider an m = 1 mode with tq = 0.25 /lis, 
w = 6.3 /im and fl p ^ = (2ir) 4 MHz. This pulse contains 



little attenuation or distortion, through the BEC, with 
the group velocity v g oc f^/IV'i | 2 - As it is input, the 
pulse induces coherent transfer of atoms from |1) to |2). 
Once the pulse is completely input (0.6 /Lts), N% nearly 

reaches Np; see Fig. E^c). Because of the finite band- 
width of the EIT transmission window there is some non- 
adiabatic loss from the BEC during the propagation, due 
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\ip2\ 2 /p, where p = |'0]P + |V'2| 2 (note the density is cylin- 
drically symmetric .23]), and cuts of the density and 
phase. The slow light propagation introduces virtually 
no phase gradient along zorr. 

For times short compared to the BEC dynamics, de- 
termined by the GPEs J5J with ^li — (typically mil- 
liseconds), any spatial mode can then be robustly stored. 
However, much longer storage can be achieved if the 
pulse parameters are chosen such that ipi and -02 remain 

(G) 

nearly stationary. In the weak probe limit ipi ~ ip\ , 
while V>2 will evolve according to an effective potential 
V2_off( r )i determined by the sum of the centrifugal po- 
tential m 2 /r 2 , the trap V(r) and the mean field potential 
f^i2|'0i| 2 ; see Fig-EIb) ■ Our choice of width w and length 
(which is governed by the pulse length To and group veloc- 
ity), are such that ip2 nearly matches the ground state of 
the potential V2-eff , which is trapping and approximately 
harmonic for an > a\2- This results in the dynamics be- 
ing greatly suppressed, as seen in the second and third 
panels of Figs.[5Jd-e). The dominant motion is a residual 
dipole sloshing in z, an unavoidable artifact of the group 
velocity and pulse length having a small dependence on 
r. For our parameters 148 ms ~ the half-period of \2-eS, 
so ip2 at this time is nearly the mirror image of ip2 at the 
switch-off time. 

After an arbitrary time, one can then switch the cou- 
pling field back on, regenerating a probe pulse accord- 
ing to fl p — — Oc n ' (V J 2 / V , i G ' ) ) Ell which is then out- 
put. Fig. |3c) shows N 2 and 7Vi oss upon a switch-on at 
148 ms. The total accumulated of losses after the pulse 
has been output is 0.33A^ in) and so A^ out) = 0.67JV"£ n) 
photons are output, nearly identical to the 63% trans- 
mission we would expect if one had propagated the LG 
mode through the BEC without stopping and storing it. 
Furthermore, the vortex phase profile has been preserved 
and little other unwanted phase gradients are introduced. 

We performed a series of numerical calculations of out- 
put pulses for the case in Fig. [21 but with switch-ons at 
various times. Since, ultimately, useful information stor- 
age requires an alphabet with several different angular 
momentum states, we also varied m. The dynamics are 
sufficiently suppressed and the variation of the storage 
efficiency T = 7Vp OUt ' ) /Np with storage time is not sig- 
nificant; see Fig. [21(f) . The variations of T can be ac- 
curately estimated by the model of [l8|. They are due 
to spatial features generated in the BECs during the dy- 
namics which get written into small temporal features 
on the revived probe fields, resulting in additional EIT 
bandwidth loss. We will report more details of these 
studies elsewhere. 

The general utility of this method for storing fields 
with a non-uniform phase profile depends not only on 
the number of photons output but also on maintaining 
the phase pattern during the storage. However, if, e.g ., 
the vortex is pinned down by an external potential |22j, 



the topological stability of the vortices in BECs makes 
the winding number m very robust, even in the presence 
of substantial unwanted phase gradients. 

We have shown how vortices can be written from 
atomic to optical fields and vice- versa, providing a unique 
nonlinear technique to create optical vortices as well as 
directly image phase singularities in BECs. This method 
should also allow the transfer of vortex configurations be- 
tween BECs which are optically connected. Furthermore, 
we considered the ability of a BEC to robustly store a vor- 
tex state input from light fields. While the case studied 
here allowed 70% efficiency, this number can be further 
improved by using larger BECs (with the loss oc 1/ ^JD\). 
We have found that an appropriate choice of parameters 
can suppress the BEC dynamics. The ultimate storage 
times should then only be limited by decoherence due to 
the small inelastic collision rates and thermal and quan- 
tum fluctuations (> 100 ms) |20j. 
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